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A three- component Fermi gas near a broad Feshbach resonance does not have a universal ground 
state due to the Thomas collapse, while it does near a narrow Feshbach resonance. We explore 
its universal phase diagram in the plane of the inverse scattering length l/akp and the resonance 
range R^kp. For a large R^kp, there exists a Lifshitz transition between superfluids with and 
without an unpaired Fermi surface as a function of l/akp- With decreasing R^kp, the Fermi surface 
coexisting with the superfluid can change smoothly from that of atoms to trimers ( "atom-trimer 
continuity"), imitating the quark-hadron continuity in a dense nuclear matter. Eventually, there 
appears a finite window in l/akp where the superfluid is completely depleted by the trimer Fermi 
gas, which gives rise to a pair of quantum critical points. The boundaries of these three quantum 
phases are determined in regions where controlled analyses are possible and are also evaluated based 
on a mean-field plus trimer model. 

PACS numbers: 03.75.Ss, 05.30.Rt, 67.85.Lm 



INTRODUCTION 

Physics of a dilute two-component Fermi gas with a 
short-range interaction becomes universal in the vicin- 
ity of a broad Feshbach resonance. In the limit of the 
vanishing potential range tq compared to the s-wave 
scattering length a and the Fermi wavelength /cp^, the 
physics is characterized solely by their dimensionless ra- 
tio l/akp. With increasing l/a/cp, a two-component 
Fermi gas exhibits a crossover from a Bardeen-Cooper- 
Schrieffer (BCS)-type superfluid to a Bose-Einstein con- 
densate (BEC) of diatomic molecules The BCS- 
BEC crossover and its variations with mass and density 
imbalances have been subject to extensive research [j42|- 

This research field can be further extended to multi- 
component Fermi gases. In particular, a three- 
component Fermi gas has attracted considerable interest 
partially because of its intriguing analogy to a matter of 
quarks with three colors [8-10], and more recently, be- 
cause of its experimental realization with ^Li atoms [lll - 
ITj . Superfluid pairing and the BCS-BEC "crossover" of 
a three-component Fermi gas have been studied theoret- 
ically 18l-[3l|, typically, in mean- field approximations. 

However, a serious problem arises when a many-body 
ground state of a three-component Fermi gas is consid- 
ered beyond the mean-field approximation: Three dis- 
tinguishable fermions can form an infinitely deep bound 
state in the zero-range limit ro ^ 0, which is known as 
the Thomas collapse [32]. Therefore, a three-component 
Fermi gas does not have a many-body ground state, or 
even if exists, it is set by a nonzero tq and hence not 
universal. This is in sharp contrast to a two-component 
Fermi gas where the Thomas collapse does not take place 
due to the Pauli exclusion principle. 

The lack of a universal ground state of a three- 
component Fermi gas can be cured in the vicinity of a 
narrow Feshbach resonance where both a and the reso- 



nance range are much larger than ro [33|. Because 
R^ sets a ground state of three fermions even in the 
zero-range limit ro ^ 0, a dilute three-component Fermi 
gas now becomes universal and its physics is completely 
characterized by the two dimensionless ratios; l/akp and 
R^kp. Recently, a two-component Fermi gas near a nar- 
row Feshbach resonance was studied experimentally us- 
ing ^Li atoms [34] and thus its extension to three com- 
ponents of fermions is quite feasible. 

The purpose of this work is to explore the univer- 
sal phase diagram of a three-component Fermi gas with 
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FIG. 1. (Color online) Minimal phase structure of a three- 
component Fermi gas with equal masses, densities, and in- 
teractions at zero temperature. There are three quantum 
phases consisting of superfluids (SF) with and without an un- 
paired Fermi surface and a Fermi gas of trimers (Trimer FG). 
The Fermi surface coexisting with the superfluid can change 
smoothly from that of atoms (FSa) to trimers (FSt)- The 
phase boundaries are asymptotically given by Eqs. (|7| and 
(|9| for the superfluid transition (solid curve) and by Eqs. (|4| 
and ()lQp for the Lifshitz transition (dashed curve). See Fig. [3] 
for an evaluation based on a mean-field plus trimer model. 
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equal masses, densities, and interactions. Our findings 
are summarized in Fig. [H There are three regions in 
the plane of l/a/cp and R^kp where quantum phases are 
easily located: a BCS superfluid with an unpaired Fermi 
surface in the weak coupling limit l/a/cp ^ — 1, a BEC 
superfluid with no unpaired Fermi surface in the strong 
coupling limit l/akp ^ 1, and a Fermi gas of deeply 
bound three fermions (trimers) in the broad resonance 
limit R^kp <^ 1. In what follows, their phase bound- 
aries will be determined in two regions where controlled 
analyses are possible: One is the narrow resonance limit 
R^kp oo with fixed l/a/cp (upper side of Fig. [T]) and 
the other is the dilute limit /cp ^ with fixed i?*/a (lower 
left and right corners). A similar approach was used in 
Ref. [35| to investigate the phase diagram of a bilayer 
Fermi gas. We also evaluate the phase boundaries based 
on a mean-field plus trimer model (see Fig. [3]) and con- 
clude that the minimal phase structure is already rich as 
shown in Fig. [TJ In particular, the Fermi surface coexist- 
ing with the superfluid can change smoothly from that of 
atoms to trimers, imitatin g th e q uark-hadron continuity 
in a dense nuclear matter |36L l37|. 



NARROW RESONANCE LIMIT 

A three-component Fermi gas near a narrow Feshbach 
resonance with SU(3) x U(l) invariance is described by 

2m 



C 



4m 



(1) 



Here ipi and represent fermionic atoms and bosonic 
molecules, respectively, and sums over repeated indices 
i, j, /c = 1, 2, 3 are assumed. The bare couplings v and g 
are related to the s-wave scattering length and the reso- 
nance range by 

V m mA 
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(2) 



A controlled analysis is possible in the narrow resonance 
limit R^ ^ oo where the Feshbach coupling g becomes 
vanishingly small [38]. In this limit, the zero-temperature 
mean- field theory yields the exact result of the systematic 
large R^kp expansion up to its leading order [39^. 

In order for all three components of fermions to have 
the same number density Ui = /c|/(67r^), three order pa- 
rameters have to be equal up to arbitrary phases. By 
introducing a gap parameter A = \/3^|(0i)|, the grand 
potential is given by 
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where Cp = /{2m) and Ep = Y^^(e^"^^70^"H^^- The 
second term is a contribution of unpaired fermions and 
the last one is that of paired fermions, both of which 
are superpositions of the original fermion components. 
For example, the unpaired fermion represented by a blue 
particle in Fig. [T] is an equal superposition; | unpaired) = 
|1) + |2) + |3). While the superfluid always exists A 7^ 0, 
the coexisting unpaired Fermi surface of atoms appears 
only for /i > and disappears for /i < 0. Therefore, a 
Lifshitz transition takes place at /i = |18j . 

By simultaneously solving the gap equation 
S^mf/SA = and the number density equation 
3ni = —dVtMF/dji^ a location of the Lifshitz transition is 
found to be 



^^^^ = "3^^^ 



(4) 



a/cF>0 



This equation determines the asymptotic behavior of the 
dashed curve in Fig. [1] toward the narrow resonance limit 
R^k^ 00. We remind that the zero-temperature mean- 
field theory (j3j) is correct only up to the leading or- 
der of the large R^k^ expansion [39|. Accordingly, the 
last term in Eq. (H]) is modified by beyond-mean-field 
corrections. The same Lifshitz transition has been ob- 
served at l/gfep = 0.633195 in the broad resonance limit 
R^kY |q|,|28|,[^, but the mean-field analysis breaks 
down here. 



DILUTE LIMIT 

A different controlled analysis is possible in the dilute 
limit /cf ^ where the problem reduces to a few-body 
problem. In vacuum, two fermions with different com- 
ponents form a bound state (dimer) with binding energy 
E2 < 0: 



-1 + v/l + 4i^,/a 



2R^ 



(5) 



a>0 



On the other hand, a binding energy of three distinguish- 
able fermions is determined by three coupled integral 
equations [4q|: 
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Two types of solutions are potentially allowed: One is 
Xi(p) = ~X2{p) with Xsip) = where Eq. (|6]) turns 
out to be equivalent to a problem of two-component 
fermions and hence no bound state solution. The other 
is x{p) = Xi{p) foi" ^ = 1^2,3 where Eq. ([6]) turns 
out to be equivalent to a problem of three identical 
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FIG. 2. (Color online) Binding energy of a ground state trimer 
E3. The dotted and dashed curves are — i?*^m|£^2| and 
— i?* y^3m|^2|/2, respectively, where E2 is the dimer binding 
energy in Eq. (|5])- 



gives 



kxT' 3 „ 
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This is a quantum critical point to separate a Fermi gas 
of trimers on its left side from the trimer Fermi gas coex- 
isting with the dimer superfluid on its right side. With 
increasing R^/a further, the Fermi surface of trimers dis- 
appears at ridimer = / {^tt'^) which gives 



-En 



a 



0.359011. 



(10) 



This is a Lifshitz transition and the system beyond this 
point is the dimer superfluid with no unpaired Fermi sur- 
face. Therefore, Eqs. (j9]) and (p!Q|) determine the asymp- 
totic behaviors of the solid and dashed curves in Fig. [TJ 
respectively, toward the dilute limit ak^^R^k^ 0. 



bosons [33I, |41|. A ground state trimer is found in the 
5-wave channel; x{p) — x(\p\)i whose binding energy 
^3 < is plotted in Fig. [2l While there is also an infi- 
nite tower of excited states due to the Efimov effect [42] , 
they are irrelevant to the present purpose to determine 
the many-body ground state. We also note that bound 
states formed with more than three fermions are unlikely 
due to the Pauli exclusion principle. 

With increasing R^/a on the side of a < 0, the trimer 
appears from the three-atom threshold at 

Es = ^ — = -0.0917249. (7) 

a 

On the left (right) side of this point, the dilute limit 
of a three-component Fermi gas reduces to a system 
of fermions (trimers). Because different components of 
fermions with equal densities can form Cooper pairs by 
an infinitesimal attraction, this is actually a quantum 
critical point to separate the atomic superfluid with an 
unpaired Fermi surface from a Fermi gas of trimers. 
Therefore, Eq. ([7|) determines the asymptotic behav- 
ior of the solid curve in Fig. [1] toward the dilute limit 
la/cpl, i^*A:F 0. 

We now turn to the side of a > but away from the 
point for Es = E2 <^ R^/a = 2.18151 where the trimer 
disappears into the atom-dimer threshold. In this region, 
the dilute limit of a three-component Fermi gas reduces 
to non-interacting dimers and trimers because their sizes 
become negligible compared to a mean interparticle dis- 
tance. Since all dimers, if exist, condense at zero tem- 
perature, the chemical potential is fixed by their binding 
energy; 2/i = E2. Then the total number density is a 
sum of contributions of condensed dimers and trimers: 



k^ 



ridh 
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With increasing R^/a, dimers appear at ndimer = which 



PHASE DIAGRAM 

Finally, we combine the above results to establish the 
universal phase diagram of a three-component Fermi gas. 
Because we developed controlled understanding on quan- 
tum phases in all available limits of l/a/cp and R^k^^ one 
emergent phase boundary has to end up as another phase 
boundary of the same type. Accordingly, the quantum 
critical lines for the superfluid transition found in Eqs. ([7]) 
and (|9]) are connected and those for the Lifshitz transi- 
tion found in Eqs. (j4]) and (p!Q|) are also connected. This 
leads to the minimal but rich phase structure shown in 
Fig.ffl 

An interesting observation is that the unpaired Fermi 
surface coexisting with the superfluid is that of atoms 
in the weak coupling l/a/cp <C —1 or dense region 
R^kp ^ 1, while it is that of trimers in the strong cou- 
pling and dilute region; l/a/cp ^ 1 and i?*/cF <C 1. 
Because there is no sharp distinction between them in 
terms of symmetries or topologies, the Fermi surfaces of 
atoms and trimers can be smoothly connected. If we 
are allowed to associate three components of fermions 
with three colors of quarks [sl-flQt , this observation corre- 
sponds to a smooth crossover from deconfined quarks to 
confined baryons with decreasing density, which is known 
as the quark- hadron continuity H, 13] • Accordingly, the 
analogous new feature in a three-component Fermi gas 
shall be termed an "atom-trimer continuity." 

Regarding low-lying excitations in the superfluid 
phase, the SU(3) x U(l) symmetry is broken down 
to SU(2) X U(l), which generates only three Nambu- 
Goldstone bosons [20]. Furthermore, in the presence 
of an unpaired Fermi surface, there exists a gapless 
fermionic excitation which has a character of atom on 
one side and that of trimer on the other side. Because 
an unpaired atom and a trimer cannot be distinguished 
by their charges in the superfluid phase, the gapless 
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FIG. 3. (Golor online) Quantum critical lines for the su- 
perfluid transition (solid curve) and the Lifshitz transition 
(dashed curve) based on a mean-field plus trimer model. The 
shaded region indicates where the smooth crossover from the 
Fermi surface of atoms (left side) to trimers (right side) takes 
place. 




FIG. 4. (Color online) Superfluid critical temperatures for 
R^kF = 0.40 (dashed), 0.381739 (dotted), and 0.36 (solid 
curve). These values of i?*/cF are above, at, and below the 
critical resonance range for opening the complete depletion of 
the superfluid by the trimer Fermi gas, which gives rise to a 
pair of quantum critical points. 



fermionic excitation can also exhibit a crossover with its 
effective mass changing smoothly from m to 3m, which 
signals the atom-trimer continuity. 

In addition to the controlled analyses, we also develop 
a model analysis to quantify the phase diagram. A guid- 
ing principle to construct the model is that it must incor- 
porate the correct asymptotic behaviors discussed above. 
The simplest possibility is just to add a contribution of 
non-interacting trimers to the mean-field grand potential 
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By simultaneously solving the gap equation 

dCluF+T/dA = and the number density equation 

Srii = —dVLMF+T/dji at A ^ 0, we find that the 
superfluid transition takes place at 



(12) 



The obtained quantum critical line is plotted in Fig. [3] 
by the solid curve, which continuously interpolates the 
correct asymptotic behaviors [Eqs. ([7|) and (|9])] in the 
weak and strong coupling limits ak^ ±0. In particular, 
we find i^^/cp = 0.288325 at l/a/cp = and its maximum 
R^kp = 0.381739 reached at l/akp = 0.314545. 

On the other hand, in the range of Fig. [3l the Lifshitz 
transition takes place at Es = 3/i where the Fermi sur- 
face of trimers disappears with increasing l/akp. The 
obtained quantum critical line is plotted by the dashed 
curve, which again yields the correct asymptotic behav- 
ior [Eq. (p!Q]) ] in the strong coupling limit ak^ +0. 
Similarly, the Fermi surface of trimers appears at £^3 = 
and that of atoms disappears at /i = 0, which correspond 
to the left and right edges of the shaded region in Fig. [3l 



respectively. While they define sharp boundaries in our 
simple treatment of trimers as point-like particles, such 
sharp boundaries need not actually exist. Indeed, a more 
elaborate model with interaction terms ~ t^ipi(j)i + h.c. 
added to the mean-field plus trimer model can incorpo- 
rate the atom-trimer continuity by hybridizing an un- 
paired atom and a trimer by condensed dimer. The 
shaded region is thus meant to indicate where the smooth 
crossover from the Fermi surface of atoms to trimers takes 
place with increasing l/a/cp. 

We also extend the mean-field plus trimer model (pTj) 
to nonzero temperature with a caveat that the mean- field 
approximation overestimates the critical temperature in 
the strong coupling region l/ak^ ^ [43|. Fig. S] shows 
the obtained superfluid critical temperature as a func- 
tion of l/a/cp for three different values of R^k^- With 
decreasing R^kp from above, the superfluid is gradually 
suppressed by the emergence of trimers. Eventually, be- 
low the critical value R^kp = 0.381739, there appears a 
finite window in 1 / a/cp where the superfluid is completely 
depleted by the trimer Fermi gas, which gives rise to a 
pair of quantum critical points. These are the same type 
of quantum critical point as that studied in the context 
of Bose- Fermi mixtures 4J, |45|] . Also the same role of 



trimer in the dilute limit was discussed in Ref. [26] with 
an explicit cutoff and hence without regard to universal- 
ity. 



SUMMARY AND DISCUSSIONS 

To summarize, a narrow Feshbach resonance sets 
the stage to investigate universal aspects of a three- 
component Fermi gas. We explored its universal phase 
diagram with equal masses, densities, and interactions 
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and found the minimal but rich phase structure shown in 
Fig. [H In particular, we observed that the Fermi surface 
coexisting with the superfluid can change smoothly from 
that of atoms to trimers, imitating the quark-hadron con- 

This new fea- 



m 



tinuity in a dense nuclear matter 
ture termed an "atom-trimer continuity" provides a novel 
analogy between ultracold atoms and quantum chromo- 
dynamics and should be investigated further. 

At absolute zero temperature, even identical fermions 
can form Cooper pairs in some partial wave channel [46] . 
An interaction between unpaired fermions in the weak 
coupling region l/a/cp <C — 1 is induced by a density 
fluctuation of the other components of fermions, while 
an interaction between unpaired trimers in the strong 
coupling region l/a/cp ^ 1 is induced by a density fluc- 
tuation of the superfluid dimers which dominates over a 
direct trimer-trimer interaction. Because both induced 
interactions cause an instability in the p-wave chan- 
nel |47j, the previously unpaired Fermi surface exhibits 
the j9-wave superfluidity regardless of whether it was of 
atoms or trimers. Therefore, the atom-trimer continu- 
ity remains intact. Here the p-wave pairing of atoms 
or trimers breaks the SU(2) x U(l) symmetry of the 
5- wave superfluid phase down to SU(2). On the other 
hand, because the s-wave superfluid phase with no un- 
paired Fermi surface retains the SU(2) x U(l) symmetry, 
the previous Lifshitz transition (dashed curve in Fig. [T]) 
becomes the p-wave superfluid transition. The trimer 
Fermi gas phase is still distinguishable from the other 
two phases because its symmetry breaking pattern is dif- 
ferent; SU(3) X U(l) SU(3) by any pairing of trimers. 

Finally, in addition to generalizing our approach to the 
case of unequal masses, densities, and interactions, the 
ideas developed here can be used to investigate universal 
aspects of a two-component Fermi gas with a large mass 
ratio > 13.6 which also lacks a universal ground state in 
the vicinity of a broad Feshbach resonance [48] . 

The author thanks J. Carlson and S. Reddy for valu- 
able discussions and a LANL Oppenheimer Fellowship 
for the support of this work. 
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